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Abstract 

We search for a light scalar particle produced in single-photon decays of the T(3S) resonance 
through the process T(3S) -^7 + ^4°, A — ► invisible. Such an object appears in Next-to-Minimal 
Supersymmetric extensions of the Standard Model, where a light CP-odd Higgs boson naturally 
couples strongly to 6-quarks. If, in addition, there exists a light, stable neutralino, decays of ^4° 
could be preferentially to an invisible final state. We search for events with a single high-energy 
photon and a large missing mass, consistent with a 2-body decay of T(3S). We find no evidence 
for such processes in a sample of 122 x 10 6 T(3S) decays collected by the BABAR collaboration 
at the PEP-II B-factory, and set 90% C.L. upper limits on the branching fraction B(T(3S) — > 
-yA°) x B(A° invisible) at (0.7 - 31) x 10~ 6 in the mass range m^o < 7.8 GeV. The results are 
preliminary. 
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1 INTRODUCTION 



The search for the origin of mass is one of the great quests in particle physics. Within the Stan- 
dard Model [T], fermion and gauge boson masses are generated by the Higgs mechanism through 
the spontaneous breaking of the electroweak symmetry. A single Standard Model Higgs boson is 
required to be heavy, with the mass constrained by direct searches to ran > 114.4 GeV [2], and by 
precision electroweak measurements to run = 1291 49 GeV [3]. 

The Standard Model and the simplest electroweak symmetry breaking scenario suffer from 
quadratic divergences in the radiative corrections to the mass parameter of the Higgs potential. 
Several theories beyond the Standard Model that regulate these divergences have been proposed. 
Supersymmetry [4] is one such model; however, in its simplest form (the Minimal Super symmetric 
Standard Model, MSSM) questions of parameter fine-tuning and "naturalness" of the Higgs mass 
scale remain. 

Theoretical efforts to solve unattractive features of MSSM often result in models that introduce 
additional Higgs fields, with one of them naturally light. For instance, the Next-to-Minimal Su- 
persymmetric Standard Model (NMSSM) [5j introduces a singlet Higgs field. A linear combination 
of this singlet state with a member of the electroweak doublet produces a CP-odd Higgs state ^4° 
whose mass is not required to be large. Direct searches typically constrain m(A°) to be below 
2mft [6] making it accessible to decays of Y resonances. An ideal place to search for such CP-odd 
Higgs would be Y — > 7A , as originally proposed by Wilczek [7]. A study of the NMSSM parameter 
space [8] predicts the branching fraction to this final state to be as high as 10 -4 . 

The decays of the light Higgs boson depend on its mass and couplings, as well as on the low- 
energy particle spectrum of the underlying theory. In certain NMSSM scenarios, particularly those 
in which the mass of the lightest supersymmetric particle (LSP) is above ra T or if ra^o < 2m T , the 
dominant decay mode of ^4° may be invisible: ^4° — > x°X°> where the neutralino x° is the LSP. The 
cleanest experimental signature of such decays is production of monochromatic single photons in 
decays Y — > 7^4°, accompanied by a significant missing energy and momentum. The photon energy 
in the Y center-of-mass (CM) H is given by 

The current best limit on the branching fraction B(Y — > jX) with X — > invisible comes from a 
measurement by the CLEO collaboration on Y(1S) [9]. The quoted limits range from 1.3 x 10 -5 
for the lightest rax (highest-energy photons) to (4-6) x 10 -4 for rax ~ 8 GeV (PDG quotes this 
result as B(Y(1S) — ► jX) < 3 x 10 for rax < 7.2 GeV [10J). There are currently no competitive 
measurements at the higher-mass Y resonances. 

In the following, we describe a search for a monochromatic peak in the missing mass distribution 
of events with a single high-energy photon. We assume that the decay width of A is negligibly 
small compared to experimental resolution, as expected [TT] for m^o sufficiently far from the mass 
of 77b [12]. Furthermore, we assume that a single A state exists in the range < m^o < 7.8 GeV; 
or if two or more states are present, they do not interfere. 



'Hereafter * denotes a CM quantity 
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2 THE BABAR DETECTOR AND DATASET 



We search for two-body transitions T(35) — > 7^4°, followed by invisible decays of ^4° in a sample of 
(121.8±1.2) x 10 6 T(3S) decays collected with the BABAR detector at the PEP-II asymmetric-energy 
e + e~ collider at the Stanford Linear Accelerator Center. The data were collected at the nominal 
CM energy E cm = 10.355 GeV. The CM frame was boosted relative to the detector approximately 
along the detector's magnetic field axis by (3 Z = 0.469. 

For characterization of the background events we also use a sample of 0.97 fb -1 collected 30 MeV 
below the T(2S) resonance, a sample of 2.6 fb" 1 collected 30 MeV below the T (35) resonance, T(45) 
decays corresponding to the integrated luminosity of 4.7 fb -1 , and 4.5 fb -1 integrated above the 
T(45) resonance. We henceforth refer to these datasets as the off-resonance sample. 

Since the BABAR detector is described in detail elsewhere [13] , only the components of the 
detector crucial to this analysis are summarized below. Charged particle tracking is provided by a 
five-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH). Photons 
and neutral pions are identified and measured using the electromagnetic calorimeter (EMC), which 
comprises 6580 thallium-doped Csl crystals. These systems are mounted inside a 1.5-T solenoidal 
superconducting magnet. The Instrumented Flux Return (IFR) forms the return yoke of the 
superconducting coil, instrumented in the central barrel region with limited streamer tubes for the 
identification of muons and the detection of clusters produced by neutral hadrons. We use the 
GEANT [14] software to simulate interactions of particles traversing the BABAR detector, taking 
into account the varying detector conditions and beam backgrounds. 

3 SINGLE PHOTON TRIGGER 

Detection of the low- multiplicity single photon events requires dedicated trigger and filter lines. 
Event processing and selection proceeds in three steps. First, the hardware-based Level-1 (LI) 
trigger accepts single-photon events if they contain at least one EMC cluster with energy above 
800 MeV (in the laboratory frame). The total LI trigger rate was typically 4-5 kHz for a combina- 
tion of 24 trigger topologies, including the single-photon line which contributed the rate of 300-400 
Hz. Second, Ll-accepted events are forwarded to a software-based Level-3 (L3) trigger, which forms 
DCH tracks and EMC clusters and makes decisions for a variety of physics signatures. Two single- 
photon L3 trigger lines were active during the data taking period. The high-energy ("HighE") line 
requires an isolated EMC cluster with CM energy E* > 2 GeV, and no tracks originating from 
the e + e~ interaction region. A subset of the data, amounting to (82.8 ± 0.8) x 10 6 T(35) decays 
and 2.6 fb" 1 collected 30 MeV below the T(35), were also processed with a low-energy ("LowE") 
single-photon trigger, which requires an EMC cluster with CM energy E* > 1 GeV, and no tracks 
originating from the e + e~ interaction region. The acceptance rate of the two single-photon L3 lines 
was up to 100 Hz. Events accepted by L3 are written to mass storage, at the rate of up to 900 Hz. 

Additional requirements are applied to the events at the reconstruction stage. We process single- 
photon events if they satisfy one of the two criteria. The "HighE" selection requires one EMC cluster 
in the event with a CM energy E* > 3 GeV and no DCH tracks with momentum p* > 1 GeV. 
The "LowE" selection requires one EMC cluster with the transverse profile consistent with an 
electromagnetic shower and a CM energy E* > 1.5 GeV, and no DCH tracks with momentum 
p* > 0.1 GeV. The two selection criteria are not mutually exclusive. 
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4 EVENT SELECTION AND YIELDS 



Due to the specifics of the online and reconstruction selections, we split the dataset into two 
broad energy ranges based on the energy of the highest-energy (in the CM frame) EMC cluster. 
The high-energy region, accepted by "HighE" L3 and reconstruction selections, corresponds to 
3.2 < E* < 5.5 GeV. The backgrounds in this region are dominated by the QED process e + e~ — ► 77, 
especially near E* = E cm /2, where the photon energy distribution for e + e~ — * 77 events peaks. 
The offline event selection is optimized to reduce this peaking background as much as possible. 

The second energy range is 2.2 < E* < 3.7 GeV, which corresponds to the "LowE" online 
selection. This region is dominated by the low-angle radiative Bhabha events e + e~ — ► e + e~7, in 
which both electron and positron miss the sensitive detector volumes. In the region 3.0 < E* < 3.7 
GeV, the tail from the e + e~ — * 77 background is significant. 

A limited number of variables are available for these very low- multiplicity event samples. We 
use the following variables to select the events of interest: 

• Photon quality: number of crystals in the EMC cluster N ciys , and transverse shower moments 
LAT and a 42 [E]. 

• Fiducial selection of the primary photons: cosine of the CM polar angle cos 9* and the 
azimuthal angle <f>* The signal photons are expected to be distributed as 1 + cos 9*, while 
the backgrounds are more strongly peaked in the forward and backward directions. 

• Extra particles in the event: we require that no charged-particle tracks are found in the SVT 
and the DCH. We also apply cuts on the energy of the second-highest photon in the event ~E\ 
(computed in CM frame), extra energy in the calorimeter S e xtra = -E-totai — E 7 , computed in 
the lab frame, and the azimuthal angle difference between the primary and the second photon 
in the event §\ — (f>*. Non-zero E extra may be present in the signal events due to machine 
backgrounds. The cut on — <p\ suppresses e + e~ — > 77 and other QED backgrounds. 

• IFR veto: we cut on the azimuthal angle difference between the primary photon and any IFR 
cluster. This variable, A<^ H , rejects the e + e~ — ► 77 events in which one of the photons is 
lost in the dead regions between the EMC crystals, but is reconstructed as an IFR cluster. 

We optimize the event selection to maximize Es /y / £ r B , where Es is the selection efficiency for the 
signal, and eb is the background efficiency. We use Monte Carlo samples generated over a broad 
range < m A o < 8 GeV of possible A masses for the signal events. We also use approximately 
10% of the available dataset as a background sample for the selection optimization. This sample is 
included in the final fit. 

In the following, we present the analysis of the data in each energy range separately. We use the 
high-energy region to measure the signal yields in the mass range < m^o < 6 GeV. We measure 
the yields in the region 6 < m^o < 7.8 GeV using the low-energy region. The overlap between the 
two regions is minimal, and the events yields are consistent in the range of m^o where the regions 
overlap. 

4.1 HIGH-ENERGY REGION 

The final selection for the energy range 3.2 < E* < 5.5 GeV is summarized in Table[H The selection 
efficiency for signal is 10-11%, depending on m^a, and is below 10~ 5 for e + e~ — > 77 events. Most 
of the signal efficiency loss occurs due to the fiducial requirements: the CM polar angle selection 
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Table 1: Selection criteria for the two regions, low and high energies. 



Variable 3.2 < E* < 5.5 GeV 2.2 < E* < 3.7 GeV 

Number of crystals in EMC cluster 20 < N CTys < 48 12 < N CTys < 36 

LAT shower shape 0.24 < LAT < 0.51 0.15 < LAT < 0.49 

ct42 shower shape 042 < 0.07 042 < 0.07 

Polar angle acceptance -0.31 < cos 9* < 0.6 -0.46 < cos 9* < 0.46 

2nd highest cluster energy (CMS) E\ < 0.2 GeV E* 2 < 0.14 GeV 

Extra photon correlation cos(02 - 4>i) > -0.95 cos(02 - 4>\) > -0.95 

Extra EMC energy (Lab) E extia < 0.1 GeV E extia < 0.22 GeV 

IFR veto cos(A^ H ) > -0.9 cos(A^ H ) > -0.95 

IFR fiducial cm(6d>Z) < 0.96 



—0.31 < cos 9* < 0.6 ensures that the second photon from e + e~ — > 77 background would hit the 
barrel regions of the EMC and the IFR, and the azimuthal requirement cos(6</>*) < 0.96 vetoes the 
dead regions between the six equal-size IFR sectors. 

We extract the yield of signal events as a function of the assumed mass m^o in the interval 
< m A o < 6 GeV by performing a series of unbinned extended maximum likelihood fits to the 
distribution of the missing mass squared 

m\ = n? r{?s) - 2E*m T{3S) (2) 

in fine steps of Am^o = 0.1 GeV. After the final selection, 955 events remain in the data sample 
in the interval —5 < m? x < 40 GeV 2 . The dominant background in this region is from e + e~ — ► 77, 
radiative Bhabha, and two-photon fusion events. The background from e + e~ — * 77 is particularly 
problematic, since its distribution peaks near m 2 x = 0. We determine the probability density 
function (PDF) of this background by selecting a high-statistics sample of on-resonance e + e~ — * 77 
events with the IFR veto removed (a total of 244,462 events). We determine the efficiency of the 
IFR veto by selecting e + e~ — > 77 events with one photon in the off-resonance sample. We find 
eiFRveto = veto /A^yy no veto = (4.5 ± 1.9) x 10~ 4 , where A" 77ve to is the number of events accepted 
with the full selection (Table [T]), and N^ inove t is the number of events accepted with the same 
selection but IFR veto removed. The uncertainty accounts for the time-dependent variation in 
EiFRveto between the different off-resonance samples. We then fix the number of expected e + e~ — > 
77 events to A r 77 = 110 ± 46. 

The background from the radiative Bhabha and two-photon fusion events is described by a 
smooth exponential function fbkg( m x ) ^ exp(cm^). The parameter c and the yield of this contin- 
uum background are left free in the fit. 

Our Monte Carlo simulations estimate that the backgrounds from the generic T(3S) decays 
or misreconstructed vector mesons produced through initial-state radiation (ISR) processes are 
negligible. The ISR processes can potentially contribute peaking backgrounds at low m x . We see 
no evidence for these extra contributions in the off-resonance sample, but also vary the peaking 
e + e~ — > 77 PDF to estimate potential systematic effects. 

The signal PDF is described by a Crystal Ball [TB] function centered around the expected 
value of m x = m 2 A0 . We determine the PDF as a function of m^o using high-statistics simulated 
samples of signal events, and we determine the uncertainty in the PDF parameters by comparing 
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Figure 1: Sample fit to the high-energy dataset (122 x 10 6 T(3S) decays). The bottom plot shows 
the data (solid points) overlaid by the full PDF curve (solid blue line), signal contribution with 
m A° = 5-2 GeV (solid red line), e + e~ — > 77 contribution (dot-dashed green line), and continuum 
background PDF (black dashed line). The top plot shows the pulls p = (data — fit)/<r(data) with 
unit error bars. 
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Figure 2: Signal yields iV s i g as a function of assumed mass m^o in the high-energy dataset. Blue 
error bars are statistical only, and the red error bars include the systematic contributions. Since 
the spacing between the points is smaller than the experimental resolution, the neighboring points 
are highly correlated. 
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the distributions of the simulated and reconstructed e + e~ — > 77 events. The resolution for signal 
events varies between a{m 2 x ) = 1.5 GeV 2 for m^o ~ to a(m 2 x ) = 0.7 GeV 2 for m A o = 8 GeV. 

A sample fit, for m A o = 5.2 GeV, is shown in Fig. [TJ This fit corresponds to the signal yield of 
^Vsig = 37 ± 15, with the statistical significance of S = \J2 ln(L max /-^o) = 2.6a, where -L max is the 
maximum value of the likelihood, and Lq is the value of the likelihood with the signal yield fixed 
to zero. No other values of m A o < 6 GeV return higher significance or likelihood. The results of 
the fits in fine steps of m A o are shown in Fig. [2j 

4.2 LOW-ENERGY REGION 

The final selection for the energy range 2.2 < E* < 3.7 GeV is summarized in Table[TJ The selection 
efficiency for signal is 20%. Most of the signal efficiency loss occurs due to the fiducial requirement 
on the CM polar angle | cos#*| < 0.46, applied to suppress the background from e + e~ — ► e + e~7, 
which rises steeply in the forward and backward directions. We restrict the photon energy range to 
avoid the region E* < 2.2 GeV where the backgrounds are excessively high and the single-photon 
trigger selection requires further investigation. 

We extract the yield of the signal events as a function of the assumed mass m A o in the range 
6 < m A o < 7.8 GeV by performing a set of unbinned extended maximum likelihood fits to the 
distribution of the missing mass squared m x in steps of Am^o = 0.025 GeV. After the final 
selection, 14,947 events remain in the data sample in the interval 30 < m? x < 62 GeV 2 . The 
dominant background in this region is from the radiative Bhabha e + e~ — > e + e~ r y events, with 
contributions from e + e~ — > 77 becoming relevant at low values of m 2 x (high photon energy). 

We parameterize the background from the radiative Bhabha events by a smooth exponential 
function /Bhabha^x) <x exp^im 2 ^ + C2m x ). The parameters c\ and C2, as well as the yield of 
e + e _ — > e + e~7 events are left free in the fit. This PDF also accounts for other radiative processes, 
such as e + e~ — ► t + t~7 and e + e~ — > fi + fi~~f. 

The e + e~ — ► 77 distribution is modeled by a sample of simulated events with an equivalent 
integrated luminosity of 22 fb — 1 . The PDF has two components: a smooth continuum with a 
turn-on of e + e~ — > 37 events starting around m 2 x = 53 GeV 2 , and a broad peak with the width 
of about 2.5 GeV 2 which accounts for a forward and backward corner of 37 phase space. The 
fraction of this peak, the fraction of the continuum e + e~ — > 3j events, and the normalization of 
the e + e~ — * 77(7) events are left free in the fit. The signal PDF is described by the same Crystal 
Ball [16j function as in the high-energy region. 

A fit for m^o = 7.275 GeV, is shown in Fig. [3l This fit corresponds to the signal yield of 
iVgjg = 119 ± 71, with statistical significance of S = y/2 ln(L max /Lo) = 1.7a. No other values in the 
range 6 < m^o < 7.8 GeV return higher significance. This fit also returns A^Bhabha = 11419 ± 441 
and A 77 = 3410 ± 451. The results of the fits in fine steps of m^o are shown in Fig. [H For each fit 
where the signal yield is allowed to vary, we also allow for the variation of the background shape 
parameters. We find that the shapes of the background PDFs are independent of the assumed m^o, 
and is also consistent with the distribution in the off-resonance sample. We observe no significant 
excess of events over the range 6 < m A o < 7.8 GeV. 

5 SYSTEMATIC UNCERTAINTIES 

The largest systematic uncertainties in the signal yield come from the estimate of the e + e~ — > 77 
peaking background yield in the high-energy region and its shape (in both energy regions). Varying 
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Figure 3: Sample fit to the low-energy dataset (83 x 10 6 T(3S) decays). The bottom plot shows 
the data (solid points) overlaid by the full PDF curve (solid blue line), signal contribution with 
m A o = 7.275 GeV (solid red line), e + e~ — ► 77 contribution (dot-dashed green line), and continuum 
background PDF (black dashed line). The top plot shows the pulls p = (data — fit)/<r(data) with 
unit error bars. 
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Figure 4: Signal yields iV s i g as a function of assumed mass m^o in the "LowE" region. Blue error 
bars are statistical only and the red error bars include the systematic contributions. Since the 
spacing between the points is smaller than the experimental resolution, the neighboring points are 
highly correlated. 
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the peaking e + e~ — * 77 background contribution by its uncertainty changes the signal yield by 
±38 events for m^o = 0, with the effect decreasing with increased m^o. The uncertainty due to 
the e + e~ — > 77 PDF is largest in the low-energy region, where it contributes up to ±70 events (for 
rrijsfi = 7.4 GeV) to the uncertainty in the signal yield. 

We determine the uncertainty in the signal PDF by comparing the data and simulated distribu- 
tions of e + e~ — ► 77 events. We correct for the differences observed, and use half of the correction as 
an estimate of the systematic uncertainty. The effect on the signal yield is generally small, except 
for the region near m^o = 7.4 GeV, where the systematic variation of the signal PDF changes the 
yield by ±64 events. Such large variation is caused by high correlation with the e + e~ — * 77 yield 
in this region. The total additive systematic uncertainty on the yield is ranges between 1 and 100 
events, depending on m^o. 

We measure the trigger and filter selection efficiency using single-photon e + e~ — > 77 and 
e + e~ — > e + e~7 events selected from a sample of unbiased randomly accepted triggers. We find 
excellent agreement with the Monte Carlo estimates of the trigger efficiency, within the systematic 
uncertainty of 0.4%. We measure the efficiency of single photon reconstruction in a large sample 
of e + e~ — > ^ + /i~7, e + e~ — ► r + r~7, and e + e~ — ► jut events, and assign a systematic uncertainty 
on the reconstruction efficiency of 2%. We assign an additional 2% systematic uncertainty on the 
single photon selection. The uncertainty on the total number of recorded T(3S) decays is estimated 
to be 1.1%. The total multiplicative error on the branching fraction is 3.1%. 

6 RESULTS AND CONCLUSIONS 

We do not observe a significant excess of events above the background in the range < m^o < 
7.8 GeV, and set upper limits on the branching fraction B(T(3S) — ► 7 A ) x B(A° — ► invisible). 
We add statistical and systematic uncertainties (which include the additive errors on the signal 
yield and multiplicative uncertainties on the signal efficiency and the number of recorded T(3S) 
decays) in quadrature. The 90% C.L. Bayesian upper limits, computed with a uniform prior and 
assuming a Gaussian likelihood function, are shown in Fig. [5] as a function of mass m^o. The limits 
range from 0.7 x 10~ 6 (at m A o = 3.0 GeV) to 31 x 10 -6 (at m^o = 7.6 GeV). These results are 
preliminary. 
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Figure 5: 90% C.L. upper limits on the branching fraction B(T(3S) — ► jA ) x B(A° — ► invisible). 
The dashed blue line shows the statistical uncertainties only, the solid red line includes the system- 
atic uncertainties. 
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